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electric field is studied. The equation for the evolution of surface electrohydrody- 
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It is well known that the tangential electric field, as opposed to the normal electric 
field, has a stabilizing effect on the free surface of a dielectric liquid. The dispersion relation 
'or the surface waves propagating along the electric field direction has the following form 

Anp p 

where k is the wave number, u is the frequency, g is the acceleration of gravity, s is the 
surface tension coefficient, p is the mass density, E is the external electric field strength, 
P = {e — lY/{e + 1), and e is the dielectric constant of a fiuid (we use the cgs electrostatic 
system of units). It is seen from this relation that, if the electric field is sufficiently strong 
{E'^ ^ l3~^y/gsp), the waves with wavelengths A = 27r/fc in the range 

sE-^p-' < A < E^Pipg)-' (1) 

will travel without dispersion, i.e., u oc k. Their propagation is described by the linear wave 
equation 

Vtt = c^Vxx, (2) 

where ri{x, t) is the elevation of the free surface, and c is the (constant) wave velosity. Its 
general solution is 

ri{x, t) = f{x + ct) + g{x — ct), 

where / and g are arbitrary functions that play the role of the Riemann invariants (they are 
constant along the corresponding characteristics). The problem of taking into account the 
system nonlinearities amounts to studying the manner in which the Riemann invariants of 
the linear wave equation will evolve under their infiuence. 

For the classical problem of long waves in shallow water, the evolution of the surface is also 
described (in the linear approximation) by the equation (121). As is known (see, for example, 
Ref. Q). the .o,U,.eanty dete.,„i„es the tendency fo. the fot,„at,o. of singularities. Sueh 
tendency can be suppressed by the linear dispersion infiuence. The asymptotic behaviour 
of the (adiabatic) Riemann invariants is then described by the Korteweg-de Vries equation. 
Below we will demonstrate that, when investigating the nonlinear waves on the surface of a 
deep fiuid layer subjected to a strong tangential electric field, we encounter the situation of a 
fundamentally different kind. The point is that, for liquids with high e values, the Riemann 
invariants of the linear wave equation, taken separately, are not distorted by nonlinearities. 
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Of course, this is not to say that the nonhnearity influence can be neglected. The nonlinearity 
is responsible for the interaction of the oppositely directed waves. 

Note that the influence of the nonlinearity on the waves propagating on the surface of a 
fluid in a tangential electric or magnetic field (these problems are similar from the math- 



ematica 
(see 



icalp 

m 



oint of view) is usually investigated in the quasi-monochromatic approximation 
and the references therein). This approach allows one to obtain immediately a 
nonlinear Schrodinger equation for the complex amplitude of travelling waves and, in so 
doing, to circumvent the difficulties associated with the necessity of analysing the nonlinear 
integro-differential equations that result from reducing the original equations of motion to 
lower dimensional equations for the surface evolution. In the present work it will be shown 
that, in the strong- field limit, such integro-differential equations can be solved and, conse- 
quently, the condition of quasi-monochromaticity can be replaced by the weaker condition 

Consider the potential motion of an ideal liquid of infinite depth in an external tangential 
electric field (the potential flow approximation holds if the characteristic wavelength A far 



exceeds the penetration depth of the vortex motion a/ J^/io, where u is the kinematic viscosity 
of the fluid [7]). We assume that, in the unperturbed state, the boundary of the liquid is a 
flat horizontal surface z = (the z axis is normal to the surface of the liquid), and the electric 
fleld vector is directed along the x axis of the Cartesian coordinate system. The function 
ri{x,t) specifles the deviation of the boundary from the plane, i.e., the liquid occupies the 
region z < r]{x,t). 

The velosity potential for an incompressible hquid $ satisfles the Laplace equation 

= 

with the following conditions at the boundary and at infinity: 

(V$)2 {e-l){VfVif'-E^) 



$ — > 0, z ^ -cx), 



where ip and (p' are the electric-field potentials in and above the liquid. The expression on 
the right-hand side of the dynamic boundary condition (nonstationary Bernoulli equation) 
is responsible for the electrostatic pressure at the interface between an ideal dielectric liquid 



and air (or vacuum) in the absence of free electric charges p. The evolution of the free 
surface is determined by the kinematic relation: 

rit = ^^~ 7]^ $2;, Z = T]. 

The electric potentials and ip' satisfy the Laplace equations 

vV = 0, VV = 0. 

Since the electric field potential and normal component of the displacement vector have to 
be continuous at the interface, we should add the following conditions at the boundary: 

ip = ip', z = ri, 

ednip = dn(p', Z = 7], 

where dn denotes the derivative along the normal to the surface z = rj. The system of 
equations is closed by the condition of the electric field uniformity at an infinite distance 
from the surface: 

ip' — * —Ex, z — s> +00, 

(f — > —Ex, z — > —00. 

The above-written equations of motion have a Hamiltonian structure, and the functions 
rj{x,t) and ip{x,t) = ^\z=r] are canonically conjugated quantities jol ll(]|: 

6H 6H 
or] dip 

where the Hamiltonian H coincides with the total energy of the system: 

H= [ ^—dxdz— [ ^ dxdz — e ( ^^Jfl dxdz. 

J 2 J Snp J 8np 

It is convenient to introduce perturbations of the electric-field potentials: 

(p' = ip>' + Ex, (p = ip + Ex. 

It can be readily shown that the perturbed potentials satisfy the Laplace equations with 
trivial conditions at infinity, 

— i> 0, z —^ +00, 
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(p ^ 0, z — oo, 
and the following conditions at the boundary: 

(p' = (p^ z = r], (3) 

e{ip, -7]^ ip^) = ip'^-r]^(f'^ - E{e -1)7]^, z = r]. (4) 

Substituting the perturbed potentials into the Hamiltonian and then integrating it by parts 
over z, we obtain: 



H 



dx, (5) 



2 8np 

where we put 4>{x,t) = (p\z='rj = 'P>'\z=rj- 

Our next problem is to express the integrand in (jj)) in terms of the functions rj and i/j. 
Let the characteristic angles of the surface slope be small, a ~ |?7a;| <^ 1 (the boundary per- 
turbation amplitude is much less than the characteristic wavelength). Then the derivatives 
of the potentials involved in (jH) and © can be expanded into a power series: 

$,|,=^ = f+9,f+V, $.1.=^ = -f+a^i^T+V, (6) 

'^x\z=ri = f+d.^f~^(f), 'Pz\z=ri = -f+d,^Hf^^(j), (7) 

= T_9,t_-V, <f%=v = f.d,Hfz'<P, (8) 
where T± = e^^^''^ , and H is the Hilbert transform: 

Hf{x) = TT-l / f{x'){x' - X)'^dx'. 



In deriving these relationships we have taken into account that fz\z=o = ='=-f^/a;L=o 
f\z=r^ = T::^{f\z=o) for harmonic functions that decay as z ^ ±00. 

Inserting (0) and (jH)) into the conditions (jH)) and and using the method of successive 
approximations, we obtain the expression for 0, correct to second order in small quantities: 
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where we have introduced the notation 7 = (e — + 1). If we substitute it together with 
© into the Hamiltonian and restrict ourselves to second- and third-order terms, we finally 
obtain: 



H 



dx — 



dx. 



The equations corresponding to this Hamiltonian are 



a 



i)t - Ht]^ = a^H{7]H7]^), + «7(Wx)x + - [liHv.y - + W^Y - + 0{a'), (9) 



rit + Hi;^ = -aH{7]Hi)^)^ - a{ri^^)^ + Oia^), (10) 
where we have passed to dimensionless variables using the transformations 

t XE-^(3-^/\A7cpY^^ t, x^Xx. 

These nonlinear nonlocal equations describe the evolution of small-amplitude surface waves. 

Applying the operator d^H to Eq. Q, differentiating Eq. (fTUj) with respect to t, and then 
subtracting the resulting expressions, we arrive at 

Vtt - v.. = oid^H [{HT^tf - l{Hv,Y - (7 - 1) (wxx + HivHv,,))] + 0{a'), (11) 

where we have used the well-known relations 2H{fHf) = {HfY ~ P and H^f = —f. 
The function ip was eliminated from the right-hand side of this equation with the help of 
the expression ip^ = Hrjt + 0{a), which follows from (jlUj) . One can see that in the linear 
approximation, i.e., in the limit a — > 0, the equation (fTT|) transforms into the linear wave 
equation of the form Q and, consequently, it describes the dispersionless wave propagation. 

Consider the influence of nonlinearity on the waves that travel along or against the 
direction of the x axis. We will seek a solution of (lllj) in the form 

ri{x,t) = f]{y,T), 

where y = x ±t and r = at. In the leading order of the expansion in the characteristic 
surface-slope angle a, we obtain the following equation for the slow evolution of the wave: 

2fjr = - l)dyHUHfiyY + ivVyy) + mkfiyy)] + 0(«), (12) 



7 



where the Hilbert transform operator acts with respect to the variable y. Unfortunately, 
this nonlinear integro-differential equation cannot be solved in the general case. One can 
assume that it describes the formation of singularities in a finite time. This process may 
be suppressed by the influence of linear dispersion which should be taken into account for a 
finite-depth fluid layer. 

Let us restrict our consideration to the degenerate case 

7(e) - 1 = 0(a), (13) 

when the equation (|T^ becomes trivial: 

fjr = 0(a). 

This means that the quadratic nonlinearities have no effect on the wave propagation, i.e., 
the nonlinearity does not distort the Riemann invariants of the linear wave equation. 

It can easily be seen that the condition ()13|) holds for dielectric liquids with sufficiently 
high values of the permittivity e. Indeed, since 7 — 1 = —2{6 — 1)^^, it can be rewritten as 

e > 0{a-^). (14) 

Let, for example, a = 0.1, i.e., the wave amplitude be an order of magnitude less than the 
characteristic wavelength, and the considered small-angle approximation is quite applicable. 
Then the condition ()14|) . i.e., the condition e > 0(10), is certainly valid for such liquids as 
ethanol {e ~ 26), nitrobenzol {e ~ 36), and deionized water {e ~ 81), which are often used 
in experiments. 

Now consider the evolution of the free surface of a liquid with high permittivity within the 
framework of the equation (lllj) describing the waves travelling in both possible directions. 
Putting 7 = 1 we, up to terms of order O(a^), obtain: 

Applying the Hilbert transform operator to both parts of this equation and introducing the 
auxiliary function v{x,t) = Hrj, we find: 

vtt - = ad^{vl - v1). (15) 

Thus, if the condition ()14|) holds true, it is possible to pass from the integro-differential 
equation |TT| to the partial differential equation. 



rjtt - Vxx = adxH 
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It is apparent that the key equation of the paper, i.e., the equation (fT^. has the following 
particular solutions: 

v = F{0, ^ = x + t, (16) 



v = G{0, C = x~t, (17) 

where F and G are arbitrary functions. These solutions correspond to the dispersionless 
propagation of nonlinear waves along (or against) the direction of the x axis. Let us study 
the interaction of these waves. 

We will seek the general solution of the equation (jlSp as a power series in the small 
parameter a: 

For leading order we obtain the linear wave equation 

% = 0, (18) 

and the equation 

% = -(^C^C)« - (^C^C)C (19) 

for order 0{a). It is meaningless to consider higher orders of the expansion since it is beyond 
the accuracy of the model The equation ()18|1 has the general (D'Alembert) solution: 

A = F{0 + G{0. 

The solution of the inhomogeneous equation ()19p is given by the following expression: 

B = F{0 + G{0 - F^G - FG^, 

where we can set F = and G = without loss of generality. An important feature of the 
described scheme for solving the equation (fTH|) is that the leading-order solution is the sum 
of the exact particular solutions (fT^ and (fT7|) . Otherwise, such procedure would give rise to 
secular terms, i.e., the terms that grow algebraically with x or t. As a result, the conditions 
of applicability of the perturbation theory would be violated for t > a^^ and |a;| > a^^. 
Thus, the solution of Eq. (|T5|l has the following form: 



v{x, t) = F{x + t) + G{x -t)- a{FG)^ + 0{a^). 
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Returning to the function 77, we find that the surface evolution is described by the expression: 

7]{x, t) = fix + t)+ g{x - t) + ad^H{Hf ■ Hg) + ©(a^), (20) 

where / = —HF and g = —HG. It is obvious from this nonhnear superposition formula that 
the interacting counter-propagating solitary waves preserve their initial shapes and phases. 

Thus, we have described the interaction of the counter-propagating waves on the surface 
of an ideal dielectric liquid in the presence of a strong tangential electrical field. It turns out 
that, even in the dispersionless case, the nonlinear evolution of the system does not result 
in the formation of singularities (e.g. wave breaking) if the condition (fT^ is satisfied. As a 
consequence, the interaction of nonlinear waves does not lead to the energy transfer towards 
small spatial scales. 

It should be noted that the results of the above investigation can be extended to the case 
of a ferrofluid in a tangential magnetic field. All one has to do is to replace the electric field 
E by the magnetic field H and the permittivity e by the magnetic permeability /x. In so 
doing the condition ()14|1 for the applicability of the solution ()2()j) turns into the inequality 
/i > 0(a;~^), which can be valid for some ferrofluids. 

An analogous approach can be applied to describe the wave propagation at the interface 
of two liquids with different permittivities ei and 62- In the limiting case when the Atwood 
number A = (pi — P2)/ (pi + P2) tends to unity (then the upper fluid has no influence on the 
boundary motion), all the above results remain valid with the only correction that e is now 
the permittivity ratio, e = ei/e2- 
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